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Abstract 

' The physical situation of the collision and subsequent interaction of plane gravitational waves 

in a Minkowski background gives rise to a well-posed characteristic initial value problem in 
which initial data are specified on the two null characteristics that define the wavefronts. 
In this paper, we analyse how the Abel transform method can be used in practice to solve 
V/~j ' this problem for the linear case in which the polarization of the two gravitational waves is 

t***" , constant and aligned. We show how the method works for some known solutions, where 

problems arise in other cases, and how the problem can always be solved in terms of an 
infinite series if the spectral functions for the initial data can be evaluated explicitly. 

■ 1 Introduction 

Many classes of explicit exact solutions are now known which model the collision and subsequent 
(3JT), interaction between plane gravitational waves with distinct wavefronts which propagate into a 



X 



Minkowski background (for a review see [|[). These solutions have invariably been obtained 
using an indirect approach which starts by finding a family of solutions of the field equations in 
the interaction region. The arbitrary parameters characterising the family are then restricted 
to specific ranges in order to satisfy the appropriate junction conditions for colliding plane 
waves. Finally, the restricted class of solutions obtained is extended backwards to determine 
the initial approaching waves that give rise to that particular family of solutions. This approach 
has enabled the general structure of this class of solutions to be understood. However, little 
progress has been made using the direct approach in which the approaching plane waves are 
specified initially and the subsequent interaction between them following their collision is then 
determined by solving the associated characteristic initial value problem. 

It is known @-@ that, in the vacuum (linear) case in which the approaching waves have 
constant and aligned polarization, the characteristic initial value problem may be solved theo- 
retically using Riemann's method. Unfortunately, this method involves integrals which cannot 
be evaluated explicitly for appropriate initial data. As an alternative, Hauser and Ernst E| have 
formulated a different method for solving the initial value problem in the linear case by making 
use of Abel transforms. Unfortunately, neither of these methods can be obviously generalised 
to the nonlinear case. However, some progress in treating this general case has recently been 
achieved. Hauser and Ernst ]|]— [Q] reformulated the corresponding nonlinear (vacuum) problem 
as a homogeneous Hilbert problem in a complex plane. Using this, they have proved that a 
solution of this problem exists, but their method has still not been used in practice to deter- 
mine an explicit solution. However, another approach to boundary or initial value problems, 
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extended to include the electromagnetic case, had previously been formulated p], |J in terms of 
the monodromy transform method. Unfortunately, problems involving singularities at the point 
of collision are encountered when applying this method to the case of colliding plane waves. 
Nevertheless, these problems can be overcome. By extending the method to include dynamical 
monodromy data, a new method for solving the initial value problem has been achieved. This 
has been outlined in and is described in detail in Jll[| . 

In the present work we concentrate on the linear case in which the polarization of the two 
initial gravitational waves is constant and aligned. We formulate the characteristic initial value 
problem for this situation in detail. We then demonstrate how the Abel transform method can 
be implemented in practice by deriving some well known solutions. Of course, to implement the 
Abel transform method to determine an explicit exact solution in the interaction region for given 
initial data, it is necessary to evaluate a number of integrals. We describe the particular diffi- 
culties that arise in attempting to implement this method for some physically interesting initial 
data. Finally, we demonstrate how the exact solution in the interaction region corresponding to 
arbitrary initial data can be constructed in terms of an infinite series involving hypergeometric 
functions. 



2 A collision of gravitational plane waves 

First, we recall that a plane gravitational wave with linear polarization and arbitrary profile 
h(u) can be described by a metric in the Brinkmann form 

ds 2 = 2dudr - dX 2 - dY 2 + h(u)(X 2 - Y 2 ) du 2 . (1) 

The wave surfaces for this wave are given by u = const. It may be noted that the profile 
function appears explicitly in the metric. This has a distinct disadvantage when the profile 
takes a distributional form as, for example, in a description of impulsive waves. For this and a 
variety of other reasons, it will be convenient to transform the metric (|l]) to the Rosen form 

ds 2 = 2dudv-P 2 dx 2 -Q 2 dy 2 , (2) 

which is always continuous, and involves two null coordinates u and v. This can be achieved 
using the transformation X = P(u)x, Y = Q(u)y, r = v + \PP'x 2 + ^QQ'y 2 , where 

P" + h(u)P = and Q" - h{u)Q = 0, (3) 

so that the only field equation is 

P" O" 

Using a natural tetrad, the only non-zero component of the Weyl tensor for the metric (|2|) is 
given by 

^ = l(^-Ol=h(n). (5) 



2\Q P 

Now consider a gravitational wave (with constant polarization) which has a distinct wave- 
front. We can label the wavefront by u = 0, so that h(u) = and the space is Minkowski for 
u < 0. The wave may then be considered to have arbitrary profile h(u) for u > 0, and be rep- 
resented by the line element (||). We can also consider a second gravitational wave of a similar 
type to propagate into the initial Minkowski space from the opposite spatial direction. This 
may be considered to have a wavefront on the null hypersurface v = using the same notation 
as in (|2|). If this wave also has constant polarization (with the polarizations of the two waves 
aligned), the region u < 0, v > can also be described by a metric of the form (|2j) but with P 
and Q being functions of v instead of u. 



In this way, we have set up initial conditions for the collision of two plane waves with constant 
aligned polarization and arbitrary profiles following the wavefronts. The space-time contains four 
distinct regions — a background Minkowski region I {u < 0, v < 0) which possesses the metric 

ds 2 = 2du dv - dx 2 - dy 2 , (6) 

two wave regions II and III (u > 0, v < and u < 0, v > respectively), and a region IV (u > 0, 
v > 0) which is subsequent to the collision of the two waves and in which they interact. The 
general structure of such space-times is well known and illustrated in figure 1. 



"Focussing" singularity 




Figure 1: The structure of colliding plane wave space-times. Region I is the background 
Minkowski space, regions II and III contain the approaching plane waves, and region IV rep- 
resents the interaction region following the collision. The singularity structure following the 
collision is known from the study of families of exact solutions and is described in Q . 

For convenience, we will denote functions and parameters associated with the region II 
with the subscript +, and functions and parameters associated with the region III with the 
subscript — . Thus, the two waves which propagate into the background have profile functions 
* 4 = h + (u) for u > 0, v < 0, and ^ = h_(v) for v > 0, u < 0. 

The problem we now address is to initially specify the initial profile functions h + (u) and 
h„(v) and then to attempt to determine the metric to which this gives rise in the interaction 
region IV. This is a well posed characteristic initial value problem for which a unique solution 
is known to exist. However, as will be seen, it is not a simple matter to construct the solution 
explicitly. 

3 Analysis of the initial data 

It is found to be convenient to continue to use two future-pointing null coordinates u and v 
throughout the space-time, with the four distinct regions being separated by two null hyper- 
surfaces u = and v = which represent the wavefronts of the approaching waves. The 
background region I (n < 0, v < 0) ahead of each wave is taken to be part of Minkowski space 
with the metric (||). 

Region II (u > 0, v < 0) contains one of the approaching plane waves. It is found to be 
convenient to represent this by the line element 

ds 2 = 2du dv - e~ u + (e v + dx 2 + e~ v + dy 2 ) (7) 

in which the metric functions U + = — log P + — log Q + , and V + = log P + — log Q + depend on u 
only. For any given initial profile function h + (u), they are constructed from the corresponding 



solutions of (|3]), and must satisfy the single equation 

2U +)UU = U +)U 2 + V +yU 2 , (8) 

which corresponds to (|J). The amplitude of the initial gravitational wave is given by 

h + {u) = -\{V + , UU - U + , U V + , U ) (9) 

In order to join this region smoothly with the metric @ across the null hypersurface u = 0, 
without introducing impulsive components in the Ricci tensor, it is necessary that U + (0) = 0, 
U+, u (0) = 0, and V + (0) = 0. It is then clear from ^) that U + and U +yU must be monotonically 
increasing functions. Hence e~ u+ is a monotonically decreasing function in this region and can 
be expanded in the form 

e -U+ = x _ ( c+n )n+ + ( u n + ) 

near the wavefront u = 0, where c + > and n + > 2 are constants. It then follows from (g) 
that, near the wavefront, V + must be such that 

V + , u 2 = 2n + (n + - l)c n + + u( n +-V + (^+- 2 )), (11) 

and the amplitude function behaves as 

h + {u) = ±y2n + (n + - 1) (n+ - 2) c n + + ' 2 u^+M' 2 + o{u^ n +l 2 ^ 2 ), (12) 

where the sign is opposite to that of V +>u . (It is assumed that all functions are sufficiently 
smooth as u — > 0.) It follows from this that the parameters n + and c + are determined uniquely 
from the initial data — specifically from the behaviour of h + (u) near the wavefront. (The case 
when n + = 2 generally contains an impulsive component on the wavefront u = 0, and needs to 
be treated a little more carefully.) 

Now, it is always possible to relabel the null hypersurfaces using an arbitrary C 2 transfor- 
mation u — > u = u(u). In practice, this freedom is almost always required in the construction 
of explicit solutions. However, it has the effect of replacing the term 2 du dv in the line element 
(|7|) by 2e~ M+<yU ^du dv . Further, in order to avoid introducing a coordinate singularity and to 
maintain continuity with the background metric in the form (^), it is necessary that du/du = 1 
on the wavefront labelled u = 0, u = 0. Using this freedom, it is possible to choose some 
particular form for e~ u+ , such as 

e~ u + = 1 - (c + u) n + . (13) 

It is significant that the power of the leading term in the expansion of the metric function 
U + cannot be changed, i.e. n + is an essential parameter which is determined by the initial 
data. Physically, this parameter determines the character of the wavefront; mathematically, it 
determines the degree of continuity of U + on the wavefront u = 0. 
Retaining complete generality at this stage it is convenient to put 

er U+ = \ +/(«). (14) 

and then to use / as a coordinate in place of u. It may be observed that this function is given 
by / = P+Q+ — \ i which is explicitly determined by the initial data. To maintain the character 
of the wavefront at u = 0, it is necessary that f(u) behaves as 

f{u) = \-{c + u) n ++o(u n +) 

near the wavefront. In fact it is always possible (although not always convenient) to use the 
above coordinate freedom to introduce a special gauge (|l3|) such that 



f{u) = \-(c + uY+. 



(15) 



The additional metric function M + , which is normally introduced with any special choice of 
gauge, is then determined by the equation 



M + f = -JJ2 + ~Hf + h)V + f 2 ^ (16) 

which replaces (||). This equation is critical in the colliding plane wave problem as it contains 
a singular term on the wavefront on which /' = 0, and this must be counterbalanced by an 
appropriate singular behaviour in the function V+. In order for M + to be continuous across the 



boundary u = 0, on which /' = 0, it is necessary that V + satisfies (11). Rewritten in terms 
of the function /, this condition can be expressed as 



lim 

/-*l/2 



(\-f)(V +f f] =2k + , (17) 



where k + = 1 — l/n + i.e. 



\ < k+ < 1. 



Region III (u < 0, v > 0), contains the other approaching plane gravitational wave with the 
same (i.e. aligned) polarization. The metric can be taken in the same form as that for region II 
but with the roles of u and v reversed. It is here convenient to put 

e- u - = l + g(v), (18) 

and to use g as a coordinate in place of v. Again, it is sometimes convenient to adopt a gauge 
such that 

g{v) = \-{c_vT-. (19) 

It may be noted that the amplitude of the two waves and the character of the wavefronts are 
determined by the constants c + , c_, n + and n_. The latter two are fixed by the initial data, but 
it is always possible to make a boost u — > cu, v —* c _1 u such that either c + or c_ can be made 
equal to unity. (Alternatively, they could be made equal.) 

The gravitational wave in this region is then specified by the functions g, V-(g) and M_(g) 
which are again related by a single equation corresponding to (16), and by the single parameter 



n_ > 2 which is contained in the form of g(v). To set consistent initial conditions (and to 
ensure that M_ is at least C° across the wavefront), it is necessary that the functions satisfy 
the condition 

\l-g)(V_ g f] =2k_, (20) 



lim 

S--1/2 



where I < k_ < 1, (or = 1 — l/n_). 

It is always possible to choose initial x,y coordinates such that the metric on the wavefronts 
prior to the collision is continuous with the background metric in the form (]6|). It is therefore 
appropriate to impose the additional initial conditions that V + (^) = 0, V_(^) = 0, M + (g) = 
and M_(±) = 0. 

Region IV (u > 0, v > 0): In the vacuum case for the collision of plane gravitational waves 
with constant aligned polarization, this region can always be described by the line element 

ds 2 = 2e- M d-u dv - e~ v \e v dx 2 + e~ v dy 2 ) (21) 

where U, V and M are functions of both u and v. The field equations immediately imply that 
e~ u satisfies a wave equation, so that we can put 

e~ u = f{u) + g{v), (22) 

where f(u) and g(v) are arbitrary functions subject to / + g > 0. However, since e~ u must be 
continuous across the boundaries u = and v = 0, f(u) and g(v) must retain the same forms as 



have already been prescribed in the initial regions. If the special gauge (|l5|) and (19) is adopted, 
it is possible to make the global definitions 

f(u) = \ - (c + u) n +G(u), g{v) = \ - (c_v) n ~Q(v) (23) 

where 0(u) is the Heaviside step function, but this prescription is not essential (and is often not 
convenient). 

It may immediately be observed that a singularity of some kind will occur in the interaction 
region on the spacelike hypersurface on which / + g = 0. It can be shown (see references in 
JIJ ) that this is generically a curvature singularity, although families of solutions with infinitely 
many parameters exist in which it corresponds to an unstable Cauchy horizon. It is referred to 
as a "focussing" singularity in figure 1. 

In complete generality, we can adopt / and g as coordinates throughout the interaction 
region, in which / < 5, 9 < \ and / + g > 0. With this, the main vacuum field equation is 

(f + 9)V fg + lV f + ±V g = 0, (24) 

which is an Euler-Poisson-Darboux equation with non-integer coefficients. For any solution 
of this equation, M can be obtained, at least in principle, by integrating the remaining field 
equations 

M f = -JJ2 + ^fTgj-h(f + 9)V f 2 (25) 

M 9 = -^ + ^)-^ + ^ 2 ( 26 ) 

In this case, the initial data is defined by the functions /, g, V+(f) and V_{g) and the solution 
of (^4|) in the interaction region V(f, g) must satisfy 

V(f,±) = V + (f), V{\,g) = V_{g). (27) 

with V(^, \) = 0. If these conditions are satisfied, the required continuity of M is assured. 
With given initial data, these conditions are sufficient to uniquely determine the solution of the 
colliding plane wave problem in the interaction region, at least up to the focussing singularity. 
We recall that, for colliding plane waves, the metric functions must be continuous across the 
wavefronts, while Vf must be unbounded on / = ^ but constrained by (17). Similarly, V g must 



be unbounded on g = ^ subject to (]20j) . However, these conditions are already included in the 
initial data functions V + (f) and V_(g). 

4 The Abel transform solution to the initial value problem 

A separable solution of equation ( |24| ) can easily be obtained in the form 

V(f,g) = A(a-f)- 1 2(a + g)- 1 2 (28) 

where A and a are arbitrary constants. Since many such functions can be superposed, we 
can replace A by A(a)y / a + \ and treat a as a "spectral" parameter and A(a) as a spectral 
amplitude function. Introducing also a second function B(a), a general class of solutions for the 
interaction region IV can be expressed in the form 

Jf vo" - /V cr + 9 J 9 Vo" + /V " - 9 

By considering the limits of this arbitrarily close to the two wavefronts as and / — ► ^ , 

it can be seen that this solution is consistent with the initial data provided 

V + (f) = /* -M=*,, V_ {g) = /* ^SLd,. (30) 
Jf V 0- - / Jg sjo-g 



As observed by Hauser and Ernst Q , these integrals are exactly those that define the Abel trans- 
form. It then follows that the spectral amplitude functions A(a) and B(a) can be determined 
explicitly from the initial data functions V + (f) and V_{g) using the inverse Abel transform 

^) = -l/» WLv, flW = _i/»*kU (31) 

n J a V J — O J a y? g - a 

When comparing (|3C| ) with (p9j), it may be observed that the solution in the interaction region 
involves the introduction of a dependence upon g in the integral involving A(a) and a dependence 
upon / in the integral involving B{a). This is an initial indication of a kind of "dynamical 
monodromy data" (as will be required for the nonlinear case and explained in [|l(J and |fi"l] ]) in 
which initial spectral data on one characteristic is modified by the presence of that on the other. 
Given initial data in the form of the specified functions /, g, V + {f) and V_(g), the derivatives 



Vl(f) and V_{g) can be computed and substituted into (31). The spectral functions A(a) and 
B{a) are then uniquely determined, and the solution throughout the interaction region is given 
by (p9|). This solves the initial value problem for the linear case, at least theoretically. In 
practice, it depends on whether or not the integrals for A(a) and B(a) in ( |31~| ) and V(f,g) in 
(|29| ) can be explicitly evaluated. 

It is possible that V + (f) and V_(g) may be expressed as a sum of distinct terms. It may then 
be observed that the Abel transform method (|29l)-(j3ll) is linear. Thus, each component of V± 
will give rise to separate components of the spectral functions A{a) and B(a). These different 
components will give rise to different components of the solution in the interaction region. Since 
the main field equation ( |24| ) is also linear, the complete solution may be obtained as a simple 
sum over all these components. 

Of course it is not necessary to construct the solution in the interaction region using a single 
integral. In cases in which the expressions change over different intervals (e.g. for sandwich 
waves, or waves with successive impulsive components), it is appropriate to divide the interaction 
region into a number of sub-regions and to determine the solution in each region successively. 
Initial data can be specified on any pair of null characteristics. Starting from any particular 
wavefront given by / = /i, the Abel transform given above can immediately be generalised to 

f fl A(a) 1 fh Vi(f) 

V + (f) = V + (h) + / da, A(a) = — / -^U= df. (32) 

This enables the Abel transform method to be used in a sequence of regions. 

In particular, the solution of the characteristic initial value problem, in which data are given 
on the two characteristics f = fi and g = g\ , is given by first determining the spectral data 

A(a) = / df, B{a) = / dg. 33) 

The solution can then be obtained by evaluating 

V{f,g) = Vifugi) + / ^^=da+ / —===d(j. (34) 

5 Illustrations of the Abel transform method 

In this section, we will illustrate how the above method can be used explicitly to construct some 
specific solutions. 



As a first example, we will derive the Khan-Penrose solution [12] which describes the collision 
of impulsive gravitational waves with aligned polarization. In this case, h + (u) = 6(u) so that 
P + = 1 — u0(u), Q + = 1 + u@(u), and the solution in region II is given by 

e -u + = i- u \ e v + = i^. (35) 

1 + u 



Thus n + = 2, and we may put 



and 



V + (f)=log 1- Jh-f -log [1 + Jh-f 



Similarly in region III, h_(v) = 5{v), and the solution is given by 



so that n_ = 2, and we can put 



9 



and 



1-v 2 



V_(g) = log 1 



l-v 
1+V 



9 



log 1 + J\-g . 



(36) 



(37) 



(38) 



We may now use the Abel transform to determine the "spectral" functions A(a) and B{a) 
that are required in the above solution of the (linear) initial value problem. In this case, using 
initial data on the null characteristic v = 0, u > 0, we obtain 



1 



vf^Jk-fa+f) 



df 



i 

7T 



: tan 




-1/2 



(39) 



By the identical calculation, using the initial data on the other null characteristic u = 0, v > 0, 
we obtain similarly that B(a) = — (i + cr) -1 / 2 . 

The solution in the interaction region may finally be obtained by evaluating 



V(f,g) 



da 



1 



9 V " + fV a ~ 9 



f ~ JV<7 + 9 
log (2a - f + g + 2^/a - fVa~T~g 



da 



log (2a + f-g + 2y/a + fV°~ 



+ log 





5 + / + 



(40) 



9, 



This is exactly the Khan-Penrose solution. To obtain the complete solution, it only remains to 
integrate the subsidiary equations (p5|) and ( p6|) to obtain the required expression for M. 

The above example is a particularly important illustration of the method for constructing 
explicit solutions for given initial data since it describes the collision of impulsive gravitational 
waves. For such waves, the metric possesses the lowest possible degree of continuity across the 
wavefronts. In this sense, the above solution represents the "worst case", and it is important to 
observe that the above Abel transform method works perfectly adequately. 

The above solution can immediately be generalised to that of Szekeres [Q] in which n± > 2 
simply by multiplying V + (f) and V_(g), and hence A(a) and B(a), by the constants ^/2k + 



and \j2k_ respectively. For this case the two terms in (40) are multiplied by \/2k + and \/2k_ 
respectively. However, for this solution, the initial profiles are given by 



h+(u) 



n 4 



5(u) 
T) (n + - 2) 



,(n+/2)-2 



2\/2 



(1 _ u n+)3-(2/n+) 



e(u) 



if n + = 2 



if n + > 2 



(41) 



with equivalent expressions for h_(v). These have the unphysical property that, apart from the 
Khan-Penrose case, the wave amplitude becomes unbounded in the initial regions as u — ► 1, 
or /--i 

The above example is a member of a family of solutions of (^4|) that can be expressed as 
products of Legendre functions. For these it is convenient to introduce new coordinates t and z 
defined by 



t = v / i-/v / i + 5+v / |-VI + / 



\-f>l\ + g->l\-9\l\ + f> (42) 

The general family of solutions is then expressed as the sum 

V = Y, (a n P n (t)P n { Z ) + q n Qn(t)Pn(z) + PnPn(t)Q n (z) + b n Q n {t)Q n {z)) (43) 
n 

where P n [x) and Q n (x) are Legendre functions of the first and second kinds respectively, and 
a n, Qn, Pn and b n are series of arbitrary constants. For example, the Szekeres solution is given 
by V = -v / 2(Vfc+ + Vk-)Qo(t)P (z) - V2(y/k+ - Vk-)P (t)Q (z), which also includes the 



Khan-Penrose solution. To satisfy the colliding plane wave conditions ( |17| ) and (20), which 
require that the derivative of V with respect to / and g be unbounded on the wavefront, it is 
necessary that at least one Legendre function of the second kind be included. 

When considering the initial data that gives rise to these solutions, the associated spectral 
functions can readily be calculated, and the solution corresponding to each component can be 
determined, at least in the simpler cases. For example, the spectral functions for the initial data 
associated with the following components are: 

V = Q (t)P (z) : A(a) = , 1 , B{a) ' 



\ + a 2^\+a 
V = Q 1 (t)P 1 (z): A(a) = =f=, B(a) 



V = Q 2 (t)P 2 (z) : A{a)=- 12 f~ 1 , B{a) '"^ " 



(44) 



4yl + a 4 \/^ +CJ 

„ t ,\ t-% / \ A/ s 20a 3 -3a nf . 20a 3 - 3a 
V = Q 3 (t)P 3 (z) : A(a) = B{a) = 

2^\ + a 2^ + a 

The complete solution corresponding to these functions can readily be evaluated. However, 
it may be observed that, for the same reasons as those outlined above, combinations of these 
solutions do not generally correspond to realistic initial data. 

Another disadvantage of the expansion (^3|) is that the / and g derivatives of all the com- 
ponents which involve Legendre functions of the second kind are unbounded on the wavefronts. 
Thus, the structure of the wavefronts are represented by an unlimited number of these terms. 
All the coefficients appear in the conditions ( |l7|) and (^C|), and so are constrained by the two 
constants n±. 



6 Some difficulties which arise in practice 

It may be noted that, for sandwich waves or sequences of impulsive waves, h + (u) is nonzero over 
a short interval only. If we denote the end of this interval by u = u%, we can consider the case 
in which h + (u) = for u > u\. Over this subsequent interval, (|3|) implies that P+(u) and 
Q + {u) must be linear in u, and it can be shown that the spectral function must take the form 



where fx = f(ui), and c is some constant. Although the initial data can be reconstructed from 
this expression using (]30|), the metric in the interaction region can only be evaluated explicitly 
using (^) for the case of the Khan-Penrose solution for which c = oo. For all other cases, 
( |29| ) gives rise to integrals of elliptic integrals of the third kind. For such situations it is most 
unlikely that the exact solution in the interaction region can be expressed using a finite number 
of expressions involving only elementary functions. 

Using (||), it may be noted that the profile function (g) can be expressed in the form 



h + (u) = -\{V +)% 



z(v +tU -u+ 



From this, it would appear to be natural to include U + as one component of V + 
spectral function for this component takes the form 



(46) 



However, the 



A(a) = - 

IT 



tan 



\ 



a 



(47) 



Although the initial data can be reconstructed from this expression, (29) again gives rise to 
integrals of elliptic integrals of the third kind which cannot be evaluated explicitly in terms of 
elementary functions. (The degenerate Ferrari-Ibahez solution [O], for which the interaction 
region is locally isometric to part of the Schwarzschild solution inside the horizon, is exactly of 
this form. It corresponds to the case in which V± is the sum of the Khan-Penrose expression 
fl3|) or (H) and U±.) 

It may be noted that (||) can also be expressed in the form 



h+(u) 



2 v +,uu 



f'v + 



2(/ + i) 



(48) 



from which it is clear that, for most initial choices of V + (u), h + (u) must be unbounded as 
/ — > — i. This feature is indeed observed for almost all the "known" solutions for colliding 
plane waves (including the degenerate Ferrari-Ibahez solution). For the initial value problem 
considered here, however, it is considered preferable that this singularity in the initial data does 
not occur. We prefer to assume that realistic expressions for V + (u) should normally take a form 
in which such a singularity is avoided. However, apart from the Khan-Penrose solution, none 
of the known families of solutions of fl24| ) for V + (u) which satisfy the condition (|l7]) also satisfy 
this regularity condition. 

These observations would seem to argue that, for collisions of plane waves with reasonably 
realistic expressions for h + (u) and h_(v), it is most unlikely that exact solutions representing 
the interaction region will be obtainable in any explicit form involving elementary functions. 



7 Series of self-similar solutions 

We may, however, recall that there exists a family of self-similar solutions of (^4|) that have been 
described in [|l4|]-[jl6| which represent waves with distinct wavefronts. These were originally 
presented in the context of different backgrounds, and in a form in which the wavefront u = 
is given by / = 0. We will now demonstrate how this approach can be applied to the colliding 
plane wave problem. 

In order to consider the wave with the wavefront u = 0, / = we replace the coordinates 
/ and g in [ 14 ] — 1 1C ] by / — A and g + \ respectively. We then introduce new coordinates r, £ in 



the interaction region defined by 



f + 9, C= 1 f { +9 , (49) 
f + 9 



In terms of the gauge (|23|), these take the form 



r - 1 - (c uYH - (c vV- C ~ 1 + {C+Ur+ - {C ~ V)n ~ (50) 

although this restriction is not necessary. Using the new coordinates, the wavefront is given by 
C = 1, and a solution of (|24| ) can be expressed in the self-similar form 

V(T,() = T k H k (() (51) 

where k is an arbitrary real (non-negative) parameter. The condition that V = on this 
wavefront is expressed by the constraint H k (l) = for k > 0. By direct substitution, it is 
found that the functions H k (Q must satisfy the linear ordinary differential equation 

(1 - C 2 )H'l + (2k - l)(H' k - k 2 H k = 0. (52) 

With the above initial conditions, these functions satisfy the recursion relations 

H k (Q = £ Hk^iOdC' so that tf£(C) = fT fc _i(0. (53) 

For integer values of k, explicit expressions for these solutions can be obtained in terms of 
elementary functions from the initial solution flo(C) = cosh -1 (. 

The solutions described above can be expressed in terms of standard hypergeometric func- 
tions F(a, b;c; z) in the form 

(/ + ,)'g, Q&) = c C "J^l +k F (U ; I + k ; L^jj (54) 
where, for integer k 

When applying the above recurrence relation for arbitrary values of k, we only require the 
recursion relation Cfe_i = — ^(k + \)c k . 

It may be noted that the condition (|17]) on the wavefront is satisfied only for the case k = 0, 
and then only with the additional multiplicative constant ao = 2^j2k + . However, in view 
of the linearity of this case, an arbitrary number of higher order terms may also be included. 
In addition, the equivalent solutions with a distinct wavefront which propagate in the opposite 
direction, and have the wavefront v = on which g = ^, can also be included. In this way, a 
general solution in the interaction region can be expressed as 

oo oo 

V(f, g) = J2 a k(f + 9?H k (±0) +$>(/ + g) k H k , (56) 

fc=0 k=0 



where ao = 2^j2k + and bo = 2y/2k_. It may be observed that these series have the advantage 
that terms with higher values of k have a higher degree of differentiability on the wavefront. This 
may be contrasted with other series representations (in terms of Legendre or Bessel functions 
etc.) in which all terms affect the differentiability properties on the wavefront. On the other 
hand, it may well not be possible to explicitly evaluate the quadratures for the remaining metric 
function M. Nevertheless, many terms can be evaluated using the methods described in |14j|— [O] 
and some qualitative properties may be deduced. 

The solution described above may be considered to be valid throughout the interaction region 
/ < i, g < i and / + g > 0. It can then be extended back to region II, according to the 



-t ^ 



Penrose construction, simply by setting g = \- The resulting expression must then coincide 
with that prescribed by the initial data. Thus, 

V(J, \) = £ a k {\ + f) k H k ( ) = V + (f). (57) 

Similarly, on the junction with region III, we require 

00 / 3 \ 

V{\,g) = Y,U\ + 9) k H h (|g )=V-(g). (58) 

The given initial data V + (f) and V_(g) can therefore be re-expressed in terms of two sequences 
of constants a k and b k . However, to address the initial value problem, we need to find a way of 
determining the two sets of constants a k and b k from the initial data functions V + (f) and V_(g) 
respectively. 

Once the coefficients a k and b k are known, the spectral functions associated with each com- 
ponent can be obtained by evaluating 

M°) = ~l£ ^==df, where V +k (f) = a k (\ + f) k H k , (59) 

and similar expressions for B k (a) and V^ k {g) involving the constant b k . It can be seen, however, 
that the integral ( |5"9j ) must take the form 

A k {a)= Pk {\-a) k (} 1 + a)- 1 l\ (60) 

where p k is a constant that is to be determined. This result has been obtained by substituting 
into the inverse transformation to reconstruct the original component. Moreover, there is no 



Pk- 



difficulty in including also an explicit ^-dependence as required in the general expression (2S). 
Explicitly, the component in the interaction region corresponding to A k (a) is 

V k (f,g) = Pk . V , / = da 

_ r(|)r(i + fe ) L t i-A 

" Pk T(l + k) yj+-g [*>*>* +k > f + g ) W 

which is exactly of the required form with 

^'-^rfr^ (62> 

The spectral functions for the complete solution of the colliding plane wave problem expressed 
in terms of these functions are therefore given by 

00 00 
A(a) = J2 Phil ~ °t{\ + °r Xj \ B(a) = £ " °)\\ + (63) 

fc=0 k=0 



where po = \/2k + and qo = \/2k_. 

Remarkably, this provides a feasible way to solve the above characteristic initial value prob- 
lem. It is not necessary to determine the two sets of constants a k and b k from the initial data 
functions V + (f) and V_(g). Rather, it is only necessary to determine the spectral functions from 



the initial data. The functions A{a)J \ + and B(a)y^ + should then be expanded as 
power series in \ — The resulting coefficients are exactly the coefficients p k and q k from which 
the coefficients a k and b k can be determined using (^) and the equivalent expression for q k . 
In this way, the general solution (|56|) in the interaction region has been constructed explicitly, 
albeit in terms of a pair of infinite series. 



-1 r\ 



8 Observations 



It has been demonstrated above that the initial value problem for colliding plane waves has 
been solved in principle for the vacuum case in which the approaching gravitational waves 
have constant and aligned polarization. The initial data are specified by two functions h + (u) 
and h_(v) which describe the profiles of the approaching waves. The first step is to solve 
the differential equations (||) for P± and Q ± . From these, the functions U + (u) and U_(v) are 
determined, and this specifies the functions f(u) and g(v). They also determine the functions 
V + (u) and V-(v), and the next step is to construct the spectral data functions A(a) and B(a) 
using (pi]). From these, the final solution for V(f,g) in the interaction region is obtained by 
evaluating the integrals in (p9|). 

The method described above is correct in principle. In practice, however, difficulties will 
occur at several stages. The greatest difficulties in the construction of exact solutions will arise 
in the explicit evaluation of the integrals for A(a) and B(a) in (|3l"|), and then particularly in 
determining the final expression for V(f,g) using (p9|). Also, although expressions for f(u) 
and g(v) are determined explicitly using the above procedure, it is not often that they can be 
inverted to construct V+(f) and V_(g). Nevertheless, it is possible to continue to use u and v as 
coordinates throughout the above calculations. 

Of course, the above procedure could be used as the basis for a numerical solution of the 
problem. Indeed, for the collision of two pairs of impulsive waves or two sandwich waves, this 
may be the only viable option. For such an approach, care needs to be taken in evaluating the 
integrals, particularly near the limits as the integrands there behave as \ f — o"| -1 / 2 , and Uf(^) 
and V_{\) are also unbounded. 

For the linear case considered above, the metric is diagonal and the solution is derived with 
the aid of the Abel transform. Unfortunately, this approach cannot be extended to the more 
general case in which the waves are not colinear. In the general case, the approaching gravita- 
tional waves either have variable or nonaligned polarization or electromagnetic wave components 
are included. In either case the main field equations are essentially nonlinear. However, the field 
equations are exactly the Ernst equations and their associated quadratures which are known 
to be integrable. Nevertheless, none of the known solution-generating techniques that are as- 
sociated with these equations are adapted to initial data on a pair of null characteristics. An 
alternative approach has therefore been developed [ 10 1 , [pj] ]. However, further work is required 



to investigate whether or not this approach can deal in practice with realistic initial data. 
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